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A very high energy 7-ray signal has been detected at the 5.5a level from 
H1426+428, an x-ray selected BL Lacertae object at a redshift of 0.129. The ob- 
ject was monitored from 1995 - 1998 with the Whipple 10m imaging atmospheric 
Cerenkov telescope as part of a general blazar survey; the results of these obser- 
vations, although not statistically significant, were consistently positive. X-ray 
observations of H1426+428 during 1999 with the BeppoSAX instrument revealed 
that the peak of its synchrotron spectrum occurs at > 100 keV, leading to the 
prediction of observable TeV emission from this object. H1426+428 was moni- 
tored extensively at the Whipple Observatory during the 1999, 2000, and 2001 
observing seasons. The strongest TeV signals were detected in 2000 and 2001. 
During 2001, an integral flux of 2.04 ± 0.35 x 10~ n cm~ 2 s" 1 above 280 GeV 
was recorded from H1426+428. The detection of H1426+428 supports the idea 
that, as also seen in Markarian 501 and 1ES2344+514, BL Lacertae objects with 
extremely high synchrotron peak frequencies produce 7-rays in the TeV range. 

Subject headings: BL Lacertae objects: individual (1ES 1426+42.8) — gamma 
rays: observations 



1. Introduction 



Blazars are the main class of Active Galactic Nuclei (AGN) detected above 100 MeV 
by the EGRET experiment on th e Compton Gamma Ray Observator y (CGRO) and by 
ground-based 7-ray observatories (IMukherjee et al.l 119971 ; IWeeked l200ll ) . They comprise a 
subclass of AGN and are characterized by a highly variable non-thermal continuum, strong 
variable optical polarization, the lack of a UV-excess (or "blue bump" ) and a core-dominated 
radio morphology. BL Lacertae (BL Lac) objects are a subclass of blazars that are notable 
for their lack of prominent emission lines. The broad-band double-humped Spectral Energy 
Distributions (SEDs) of BL Lacs identified in x-ray surveys differ significantly from the SEDs 
of those identified in radio surveys. This led to the sub-classification of BL Lacs into High- 
frequency-peaked BL Lacs (HBLs) and Low -frequency-peaked BL Lacs (LBLs) based on the 
ratio of their x-ray to radio flux densities (jPadovani fc Giommilll995l ). The first "hump", 
generally assumed to be the peak of the synchrotron emission (in a vF v representation), is 
in the IR-optical for LBLs and in the EUV-soft x-ray band for HBLs. BL Lac objects make 
up a significant fraction of the 70 blazars in the 3rd EGRET Catalog (IHartman et al.lll999l ) 
and most of them are classified as LBLs. It has been shown that there is not a shar p division 
between these two classes of objects (see, e.g., iFossati et al.lll998l ; lGhisellinilll999l ). 



Recent studies with ground-based 7-ray telescopes have produced evidence for TeV 7- 
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ray e mission from seven BL Lac objects, five of which are classified as HBL and two as 
LBL (jWeekesll200ll ). The most prominent of these are Markarian 421 (Mrk 421), which has 
been detected by five ground based imaging atmospheric Cerenkov 7-ray observatories, and 
Markarian 501 (Mrk 501), which has been detected by six such observatories. The emis- 
sion from both of these objects can be explained by Compton-synchrotron models although 
detailed modeling is still fraught with many uncertainties. Both objects are characterized 
by rapid variability on time-scales from hours to months. In the TeV energy range, the 
energy spectrum of Mrk 501 and, more recently, that of Mrk 421, have been shown to ex- 



hibit absorption-like features (IKrennrich et al.lll999l ; lAharonian et al.lll999t iKrennrich et al. 



200ll ) . The temporal and spectral properties of the other TeV BL Lacs are less well defined. 



Since 1992, the Whipple Gamma Ray collaboration, using the 10m imaging atmo- 
spheric Cerenkov telescope on Mt. Hopkins, has been searching for TeV 7-ray emission 
from AGN. Initially the search was concentrated on blaz ars detected by E GRET at any red- 
shift; these observations l ed to the detection o f Mrk 421 (IPunch et al.lll992l ) and upper limits 
on some 30 other blazars (IKerrick et al.lll995l ). More r ecently, the search has concentrated on 
nearby BL Lacs leadi ng to the detection of Mrk 501 f Quinn et aDll996l ) and 1ES 2344+514 
( iCatanese et al.lll998l ). Between 1995 and 1998, the survey included 24 objects (17 HBLs, 7 
LBLs) ranging in redshift from 0.046 to 0.44; the results for these objects will be published 
shortly (D. Horan et al., in preparation). Although none of these observations resulted in a 
detection, the observations of H1426+428 yielded the highest consistently positive statistical 
significances. 

In this paper, evidence is presented for the detection of very high energy (VHE) 7-rays 
from H1426+428. The H1426+428 observations are divided into two categories - those taken 
as part of the general blazar survey between 1995 and 1998, and the subsequent concentrated 
observations carried out between 1999 and 2001. A preliminary energy spectrum is derived 
and the implications of the detection of H1426+428 at these energies are discussed. 



2. H1426+428 at Other Wavelengths 



H1426+428 was discovered in the 2-6 keV band by HEA P- 1 (IWood et al.lll984f ) and was 
classified as a BL Lac object in 1989 (IRemillard et al.lll989l ). It has a cosmological redshift 
of 0.129 and is located in the constellation of Bootes (a J20 oo = ^ h 28 m 32 s . 7, S. m0 o = +42° 
40' 20"). It is an optically faint object (my = 16.9) and is believed to be at the center of an 
elliptical galaxy. H1426+428 is bright in the x-ray band with a 2-6 keV luminosity of ~ 10 44 
erg s _1 which is typical of the BL Lacs found with HEAO-1. Both the flux of H1426+428 in 
the 2-10 keV band and its spectral index above 2 keV have been found to change over time 
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( jCostamante et al.ll200ll ). 



Three recent BL Lac surveys, DXRB S (jPerlman et al.lll998l ). RGB (ILaurent-Muehleisen et al. 



19991 ). and REX (jCaccianiga et al.lll999l ). have shown that, despite the conventional subdivi- 
sion into HBLs and LBLs, BL Lacs actually form a continuous class with respect to the peak 
of the synchrotron emission, which smoothly ranges between IR and soft x-ray frequencies 
and up to the 2-1 keV band for source s like Mrk 421. Beppo-SAX observati ons of Mrk 501 
Jpian et al.lll998h and 1ES 2344+514 jGiommi. Padovani. & Permian! 120001 1 revealed that, 
at least in a flaring state, the first peak of the SED can reach even higher energies, at or 
above 100 keV. 

In 1998-1999, Beppo-SAX performed an observing campaign with the aim of finding 
and studying other sources as "extreme" as Mrk 501 is in its flaring state. The candidates 
for the Beppo-SAX survey were selected from the Einstein Slew Survey and the RASSBSC 
catalogs. These Beppo-SAX observations (Costamante et al. 2000a, 2000b, 2001) revealed 
four new "extreme" HBLs, selected because they have high synchrotron peak frequencies 
and are therefore possible TeV emitters. These four candidates for TeV emission were: 1ES 
0120+340, PKS 0548-322, 1ES 1426+42.8 (i.e. H1426+428) and H2356-309. 

The spectra for three of these objects (1ES 0120+340, PKS 0548-322, and H2356-309) 
were well fitted by a convex broken power law with a break energy, and hence the peak of the 
synchrotron emission, occurring at about 1.4 keV for 1ES0120+340, 4.4 keV for PKS0548- 
322 and 1.8 keV for H2356-309. In the case of H1426+428 however, no evidence for a spectral 
break up to 100 keV was found. Instead, its spectrum was well fitted by a single power-law, 
with a flat spectral index of 0.92 up to 100 keV, thus constraining the peak of the synchrotron 
emission to lie near or above this value. T he best fit of a pure hom ogeneous synchrotron 
self-Compton (SSC) model for H1426+428 (ICostamante et al.ll2000al ) predicted detectable 
7-ray emission at TeV energies. At the time of the Beppo-SAX observation the observed x- 
ray flux in the 2-10 keV band was at one of t he lowest levels ever reco rded from H1426+428, 
indicating that it was not in a flaring state (ICostamante et al.l 120011 ). This implies that, in 
the event of a flare, the synchrotron peak could shift to even higher values as was observed 
for both Mrk 501 and 1ES 2344+514. It also indicates that highly relativistic electrons are 
present which makes H1426+428 the most promising candidate for TeV emission from this 
survey. 
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3. Observing Technique 



The observations reported in this paper were taken with the 10m reflector at the Whipple 
Observatory on Mount Hopkins in southern Arizona (elevation 2.3 km) using the atmospheric 
Cerenkov imaging technique. During the course of the observations p resented here, many 



changes were implemented on the 10m telescope. Th e imaging camera (jCawley et al.l Il990l ) 



1. 



was upgraded a number of times (iFinley et al J 12000), the number of photo-multiplier tubes 
(PMTs) being increased w i th ea ch iteration, resulting in the current high resolution, 490 
pixel camera (IFinlev et al.l 120011). The triggering electronics were upgraded to a Pattern 
Selection Trigger (IBradbury et al.lll999l ) and light cones were installed in front of the PMTs 
to increase their light collection efficiency. The unprotected, anodized, front- aluminized 
mirrors on the 10m were recoated during the timespan of the observations. Because of 
changes in PMT configuration, triggering, mirror reflectivity, and light collection efficiency, 
the energy response of the camera for 7-ray detection varied. Since various changes were 
made at the beginning of each observing season (typically between June and September), it is 
convenient to consider each observing season separately and to characterize it by the observed 
response from the assumed standard candle, the Crab Nebula. The different configurations 
of the camera and the resultant peak response energies are summarized in Table [TJ 
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Table 1. Camera configurations from 1995 - 2001. 



Period 


No. of 


Spacing 


FOV a 


Light 


Trigger 


P b 

E->peak 






Pixels 


(°) 


(°) 


Cones 




(GeV) 


1995/01 - 


1996/12 


109 


0.259 


3.0 


yes 


majority 


300 


1997/01 - 


1997/06 


151 


0.259 


3.3 


yes 


majority 


350 


1997/09 - 


1998/12 


331 


0.24 


4.8 


no 


majority 


500 


1999/03 - 


1999/06 


331 


0.24 


4.8 


yes 


pattern 


500 


1999/09 - 


2000/07 


490 


0.12 c 


3.8 d 


yes 


pattern 


430 


2000/10 - 


2001/06 


490 


0.12 c 


3.8 d 


yes 


pattern 


390 



a Field of view (FOV). 

b The peak response energy (E pea k); this is the energy at which the col- 
lection area folded with an E -2 ' 5 spectrum reaches a maximum. Thus, it is 
the energy at which the camera is most efficient at detecting 7-rays. These 
values are subject to a ~ 20% uncertainty. Although E pea k has increased 
somewhat over time, this does not mean that the camera is now poorer at 
detecting the lower energy 7-rays. In fact, the collection area of the tele- 
scope at 300 GeV in the 1999 - 2000 observing season, was greater than the 
collection area at this energy in the 1995 - 1996 observing season. 

c The spacing between the outer tubes is 0.24°. 

d The outer rings of tubes were not used in this analysis, hence the field 
of view here is effectively 2.6°. 
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The Cerenkov light im ages from each shower are analyzed off-line using a moment anal- 



ysis ( [Reynolds et al.lll993l ). The derived image parameters are used to distinguish candidate 
7-rays from the large background of cosmic-rays. These parameters include length, width, 
distance (from the optic axis), and orientation (alpha). In addition, the two highest recorded 
signals in individual pixels (sigl and sig2) are noted as well as the total light in the image 
(size). Monte Carlo simulations were used to determine the approximate limits on the pa- 
rameters to be used for the identification of candidate 7-ray images; these limits were then 
optimized on contemporaneous observations of the Crab Nebula. The results of an analysis 
are graphically presented as a histogram of the alpha parameter - usually referred to as an 
alpha plot. In such histograms, the alpha values for all the events that passed all 7-ray 
selection criteria except for the alpha cut, are plotted. For a 7-ray source, an excess would 
be expected at low values of alpha. 

Two modes of observation were used: on/off and tracking. In on/off mode, a 28 
sidereal minute run is taken with the candidate 7-ray source at the center of the field of view 
- the ON run. The OFF run, a 'control' run, is also taken for 28 sidereal minutes through 
the same range of azimuth and elevation, thus enabling this region of the atmosphere to be 
characterised in the absence of the 7-ray candidate. In the tracking mode, the object is 
continuously tracked for 28 sidereal minutes with no control run being taken. 

Even in the absence of a 7-ray source, a certain percentage of events recorded will pass all 
of the 7-ray selection criteria. This background level of events, which depends on a number 
of factors including sky brightness, elevation, and weather, needs to be established in order 
to calculate the statistical significance of any apparent excess of 7-rays detected. There are 
different methods to estimate this background depending on which mode of observation was 
used. 

For data taken in on/off mode, the OFF scan provides an estimate of the number of 
7-ray like events that would be recorded from this range of azimuth and elevation, in the 
absence of the candidate 7-ray source. In this data acquisition mode, differences in night-sky 
brightness between the ON and OFF regions of the sky can introduce a bias when the data are 
analyzed. In order to reduce this bias, the standard deviations of the night-sky background 
in each PMT are compared for the ON and the OFF runs. Gaussian noise is then added in 
quadrature to the signal from whichever region of the sky (ON or OFF) is the darker, so as 
to match the standard deviations for the tube in the ON and OFF runs. This technique is 



known as software padding and is described in ICawley! (119931 ) . 



Unlike observations taken in the on/off mode, scans taken in the tracking mode do 
not have independent control data. A background estimate is essential in order to predict 
what number of the background events, passing all cuts, would have been detected during the 
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scan in the absence of the candidate 7-ray source. Since most of the H1426+428 observations 
were taken in the tracking mode, two different methods of background estimation were 
used. Each of these methods is outlined below. 

In the first method, which is the s tandard analysis method used for data taken in the 



tracking mode (ICatanese et al.lll998l ). the alpha plot was characterized when there was no 
7-ray source at the center of the field of view. This was done by analyzing 'darkfield data' 
which consisted of data taken in the off mode, and of observations of objects found not to 
be sources of 7-rays. A large database of these scans were analyzed, giving the shape that 
an alpha plot would be expected to have when no 7-rays are present. 

Most 7-rays from an object at the center of the field of view will have small values 
of the alpha parameter. Hence, the alpha distribution beyond 20° can be assumed to be 
independent of the 7-ray source, and thus representative of the background level of 7-ray-like 
events in the field of view. Using the darkfield data, a ratio was calculated to scale the 
number of events between 20° and 65° to the number that pass the alpha cut. This ratio, 
the 'tracking ratio' (p), was then used to scale the 20° to 65° region of the alpha plot for the 
tracking scan, to estimate the background level of events passing all cuts. Its value and 
statistical uncertainty (Ap) are given by: 



1 a N alpha / N alpha ^ 1 alpha /-, \ 

P±Af)= N ± \W + /V^ (1) 

1 ^ control y JV control ly control 

where N a i p ha is defined as the number of events in the darkfield data that pass all the 7- 
ray selection criteria including the alpha cut, while N contro i is the number of such events 
with alpha between 20° and 65°. The tracking ratio was calculated using the darkfield data 
available for each observing year over the same range of elevation angles as the source ob- 
servations. In addition, the tracking ratio was checked using the OFF data on H1426+428 
from on/off runs (where available); within statistics these tracking ratios were consistent 
with the standard tracking ratio derived from the full yearly database. To check for system- 
atic effects associated with the determination of the global tracking ratio, the darkfield data 
were subdivided on the basis of the region of the sky which they were from. Each of these 
dark-sky regions was checked individually and the tracking ratio determined and compared 
to the global tracking ratio. It was found that the global tracking ratio and these individual 
tracking ratios were consistent within statistical uncertainty. 

In the alternative method for estimating the background, each H1426+428 tracking 
scan was assigned a suitable OFF scan as its background. The off database comprised all 
OFF scans taken during the observing season including H1426+428 OFF runs. In order to 
establish the most suitable OFF scan for each H1426+428 tracking scan, each off run was 
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characterized by parameters pertinent to the conditions during the scan. These consisted of 
the elevation, the number of PMTs switched off and the mean night-sky background during 
that scan, along with the throughput factor recorded on the night of th e scan and the date 



on which the scan was taken. The throughput factor ({Holder et al.ll200ll ) is a measure of sky 
clarity based on the spectrum of the total amount of light produced by background cosmic 
rays. The same information was assembled for each H1426+428 tracking scan. Then, 
for each H1426+428 tracking scan, the OFF scan whose conditions most closely matched 
it was deemed the most suitable background estimate, and was used as the background 
for that H1426+428 tracking run. In cases where there was not a suitable off run, the 
H1426+428 tracking run was omitted. With each H1426+428 tracking scan then having 
an associated off run, the data were analyzed as if they were taken in the ON / off mode, 
and hence, software padding was applied. The sky quality changes from night to night, 
causing variations of ~ 10% in the raw trigger rate, for data taken on the same source at the 
same elevation. Therefore, even after selecting off scans with run conditions which match 
those of the ON scans, it is still necessary to normalize the total number of events ON and 
off . We calculate the normalizing factor from the ratio of the number of events in the 30° 
- 90° control region of the ON and OFF alpha plots. 

It was found that both methods of background estimation were consistent, and indi- 
cated excesses of similar significance, thus suggesting that the 7-ray excess can be reliably 
determined for tracking scans. 



4. Database 

For the purpose of data analysis, two different strategies were employed: (a) All of the 
data taken in the ON and the tracking modes were combined and subjected to a tracking 
analysis, (b) All of the data taken in the ON mode, along with their corresponding off data, 
were analyzed together using the pairs analysis. The 7-ray selection criteria, including the 
alpha bound, were re-optimized each year using data from the Crab Nebula. The optimum 
value for the alpha bound was found empirically to increase from 10° to 15° when the smaller, 
high resolution 490 pixel camera was installed in the summer of 1999. 



4.1. 1995 - 1998 Observations 



In June- July, 1995, H1426+428 was observed for 3.5 hours in the tracking mode with 
the 109 pixel camera (spacing 0.259°, 3.0° field of view and standard trigger). In February- 
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June, 1997, another 13.2 hours were used to observe H1426+428 in the tracking mode; 
the camera configuration was as before but the number of pixels was increased to 151. In 
April, 1998, 0.87 hours of H1426+428 data were taken; the 331 pixel camera was installed 
at this time with the standard trigger. 

4.2. 1999 Observations 

During the 1999 observing season, a total of 7 on/off pairs and 51 tracking runs 
were taken with the 2-fold Pattern Selection Trigger. Three tracking runs were excluded 
from further analysis based on fluctuations in the raw rate and inferior weather conditions. 
This left a total of 24.35 hours (7 ON runs and 48 tracking runs) of data available for 
analysis. These observations were initiated as part of the general BL Lac survey (D. Horan 
et al., in preparation) with extra data being taken in 1999 because of an initial indication 
of a signal from H1426+428. A tracking ratio of 0.232 ± 0.005 was determined for alpha < 
10°. 



4.3. 2000 Observations 

During the 2000 observing season, 33 ON / off pairs and 35 tracking runs were taken 
on H1426+428. Eleven ON source runs were excluded from further analysis based on fluctu- 
ations in the raw rate and inferior weather conditions, leaving a total of 26.46 hours of data 
available for tracking analysis (57 tracking/on) and 13.86 hours of data in the on/off 
mode (30 on/off pairs). A tracking ratio of 0.312 ± 0.002 for alpha < 15° was determined. 
H1426+428 OFF data were analyzed to ensure that this tracking ratio was consistent with 
the tracking ratio derived for these data. 

4.4. 2001 Observations 

A total of 42.24 hours were spent on the source during 2001; this comprised 39 on/off 
pairs and 59 tracking runs. After excluding data taken at low elevation and in unsatisfac- 
tory weather conditions, 38.10 hours of data remained (34 on/off pairs and 53 tracking 
runs). The tracking ratio for this observing season, again for alpha < 15°, was calculated to 
be 0.323 ± 0.002. As with the 2000 data, checks were performed to ensure that this ratio 
was appropriate for the H1426+428 data. A summary of the data taken and analyzed in the 
2000 and 2001 observing seasons is give in Table [2j 
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5. Results 

5.1. 1995 - 1998 Observations 

The observations of H1426+428 between 1995 and 1998 were analyzed as part of a 
general BL Lac survey. The results are summarized in Table [31 Although not statistically 
significant, H1426+428 was one of only two objects, out of the 24 objects observed in this 
survey, to show consistently positive results. 
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Table 2: Summary of the number of data scans taken and analyzed in 2000 and 2001. 



Year 


Data Recorded 


Data Analyzed 






TRACKING ON 


Standard Alternative 


PAIRS 


2000 


35 (27) 33 (30) 


57 57 


30 


2001 


59 (53) 39 (34) 


87 71 


34 



Note. — The mode that the data were recorded in are given in the first two columns. The numbers 
in parentheses here refer to the number of usable data scans recorded. The final three columns give the 
number of scans analyzed in each of the three data analysis modes: standard tracking analysis, alternative 
TRACKING analysis or pairs. 
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Table 3. Summary of the Observations and Results for H1426+428 between 1995 and 

1999. 





Exp. 


Total 


Max. a 


Max. a 


Flux c 




Period 


(hrs) 


a 


Month a 


Night b 


(cm -2 s _1 


) 


1995/06 - 1995/07 


3.48 


2.1 


2.1 


2.1 


<2.2 x 10" 


n 


1997/02 - 1997/06 


13.16 


1.7 


2.2 


1.6 


<1.6 x 10" 


n 


1998/04 


0.87 


1.7 


1.7 


2.0 


<6.7 x 10- 


n 


1999/03 - 1999/06 


24.35 


0.9 


1.6 


2.1 


<6.7x 10- 


li 



a The maximum statistical significance of the signal recorded from 
H1426+428 when the data are grouped by the month during which 
they were recorded. 

b The maximum statistical significance of the signal recorded from 
H1426+428 when the data are grouped by the night on which they 
were recorded. 

c The integral fluxes are quoted at the peak response energy for the 
observation period, as given in Table [TJ 
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5.2. 1999 Observations 

A total of 24.35 hours of data were taken on H1426+428 during 1999. The net excess 
of +0.9cr from the combined 1999 observations was not statistically significant. 

5.3. 2000 Observations 

A tracking analysis of all the on-source data taken on H1426+428 during 2000, us- 
ing the standard tracking ratio calculated for that observing season, gave a +4.2<r excess 
corresponding to a 7-ray rate of 0.24 ± 0.06 7 min -1 . When these data were analyzed in- 
dependently using the alternative method of background estimation, the excess was at the 
+3.1cr level with a 7-ray rate of 0.21 ± 0.07 7 min -1 . The alpha plot for these data is 
shown in Figure [T] (left panel) along with the alpha plot from the matched OFF data. The 
background shown on the figure (dashed line) consists of the OFF data which were chosen 
to match the characteristics of the H1426+428 tracking data. These data are scaled to 
match the H1426+428 data in the 30° - 90° region of the alpha plot. The net significance for 
the 30 on/off pairs taken in 2000 was +1.2cr, with a 7-ray rate of 0.10 ± 0.09 7 min" 1 . The 
alpha plot for these data is shown in Figured] (right panel). On the night of May 30, 2000 
(MJD 51694), the first tracking scan on H1426+428 gave a significance of +3. Oct, using 
the standard tracking ratio method of background estimation. This prompted observers to 
take four more tracking scans which led to the detection of a signal at the +3.7cr for this 
night. The 7-ray rate for these 5 scans (2.15 hours) was 0.67 ± 0.18 7 min -1 . 
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Fig. 1.— Left: Alpha plot for 26.49 hours of the data taken on H1426+428 during 2000. The 
total significance for these 57 scans is +3. la when analyzed using the alternative method for 
estimating the background. The matched off data used to estimate the background level 
are also shown (dashed line). These are normalized to the H1426+428 data between alpha 
values of 30° and 90°. Right: Alpha plot for the 30 on/off pairs taken in 2000. The net 
significance for these data is +1.2a. 
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5.4. 2001 Observations 

A tracking analysis of all of the on-source data taken on H1426+428 during 2001, 
using the tracking ratio calculated for this observing season, gave an excess of +5.4cr corre- 
sponding to a 7-ray rate of 0.36 ± 0.07 7 min -1 . An independent analysis of 31.13 hours 
of these data, using the alternative method of background estimation, resulted in a +5.5cr 
excess with a 7-ray rate of 0.44 ± 0.08 7 min -1 . The alpha plot for these data is shown in 
Figure [2] (left panel) along with the alpha plot from the matched OFF data. As before, these 
data are scaled to the H1426+428 data between alpha values of 30° and 90°. 
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Fig. 2.— Left: Alpha plot for 31.13 hours of the data taken on H1426+428 during 2001. 
The total significance for these 71 scans, when analyzed using the alternative method of 
background estimation, is +5. 5a. The matched OFF data used to estimate the background 
level are also shown (dashed line). These are normalized to the H1426+428 data between 
alpha values of 30° and 90°. Right: Alpha plot for the 15.53 hours of data taken in the 
on/off mode on H1426+428 during 2001. The total significance for these 34 on/off pairs 
is +2.0a. 
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The significance from the 34 on/off pairs taken during the 2001 observing season is 
+2. Oct, with a 7-ray rate of 0.22 ± 0.11 7 min -1 . The alpha plot for these data is shown in 
Figure [2] (right panel). A summary of the results of the observations made on H1426+428 
during 2000 and 2001, using the alternative method of background estimation, is presented 
in Table H 
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Table 4. Summary of the results of H1426+428 observations during 2000 and 2001. 



Period 


Exp. (hrs) 


Total a 


Max. a Month a 


2000/02 - 2000/06 


26.37 


3.1 


3.4 


2001/01 - 2001/06 


31.12 


5.5 


5.0 



a The maximum statistical significance of the signal recorded 
from H1426+428 when the data are grouped by the month 
during which they were recorded. 
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When both the differences in exposure time and methods of background estimation are 
taken into account, the results from the on/off pairs and the tracking data taken on 
H1426+428 during each observing season are consistent, at the 2a level, with each other. 



5.5. Comparison of the Gamma-Ray and X-Ray Flux from H1426+428 

The 7-ray rates for the 2000 and 2001 H1426+428 data were compared with the x-ray 
flux fro m the All Sky Monit or (ASM) instrument on board the Rossi X-ray Timing Explorer 



(RXTE; iLevine et al.l 119961 ) . A correlation was sought between the nightly 7-ray rates and 
the "one-day average" x-ray data points. Only nights on which there were both x-ray and 
7-ray data were considered when performing the correlation. This left 37 nights for analysis 
during 2000, and 39 nights during 2001. 

Since 7-ray observations can only be taken on moonless nights, there are a few (~ 6) 
nights around the time of full moon each month, when no observations can be made. The 
periods during which the 7-ray data are taken are referred to as 'darkruns'. Correlations 
between the x-ray and 7-ray rates were sought both over the entire observing season, and 
during darkruns with more than two nights on which both x-ray and 7-ray data were taken. 
No evidence for significant nightly correlation was found for either observing season. The 
x-ray and 7-ray rate curves are shown in Figure [3J The rates plotted here are the average 
monthly rates for the 2000 and 2001 data. The 7-ray rates shown were calculated using the 
alternative method of background estimation. There is some evidence for an anti-correlation 
between these average monthly rates, especially in the 2001 data. 
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Fig. 3. — Top: The mean x-ray rate for H1426+428 for each month during 2000 and 2001 
from the ASM on RXTE. The data are only plotted for months during which 7-ray data 
were also taken. Bottom: The mean 7-ray rate from H1426+428 for each month calculated 
using the alternative method for background estimation. This rate is plotted in units of the 
7-ray rate from the Crab for that year. The detected Crab rate was 2.45 7 min -1 during 
2000 and 3.36 7 min" 1 during 2001. 
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6. Spectral Characteristics of H1426+428 

The TeV flux from H1426+428 is weak and near the threshold of sensitivity of the 10m 
7-ray telescope. The photon flux i s so small that i t is imp ractical to apply standard spectral 



analysis techniques, for example, iMohanty et al.l (119981 ). to these data. In fact, only the 



differential flux from the source for the energy at which the telescope is most efficient at 
detecting 7-rays, can be evaluated with reasonable accuracy. Such a flux estimate accounts 
mostly for the photons with energies corresponding to the peak in the differential detection 
rate of the telescope. This energy, E p , depends on the spectral index of the source in 
question. For observations of H1426+428 with the Whipple telescope in 2001, E p was found 
to be between 280 - 360 GeV. The rate of change of the number of excess events, presumably 
photons from the source, with the total amount of light in the Cerenkov image, the size, is 
directly related to the spectral index of the source. This relationship is established using 
Monte Carlo simulations; the spectral index of the source in the vicinity of E p is then 
evaluated. This in turn allows estimates of E p for this spectral index, and the differential 
and integral fluxes to be improved by reducing the uncertainty due to the unknown spectral 
index. 

Figure H] shows the excess events detected from the direction of H1426+428 as a function 
of integrated Cerenkov light in the shower image. For comparison, this plot also shows 
the excess events detected from the Crab Nebula after a 4.1 hour exposure, which was 
chosen to produce a similar maximal excess of ~ 10 3 events as was found during the 32.5 
hours of observation of H1426+428 during 2001. The lowest size cut of 366 digital counts 
is the result of signal to noise optimization of the H1426+428 data, producing an overall 
significance close to +6a. The highest size cut, 840 digital counts, is limited by the rapidly 
declining significance to the level of +3. 5a. To a first approximation, the dependence of 
the excess events on the logarithm of the size cut applied, is a linear function whose slope 
is monotonically related to the spectral index of the source. This plot indicates that the 
H1426+428 spectrum is softer than the well-studied spectrum of the Crab Nebula. 

Monte Carlo simulations were used to relate the observed parameters of linear fits to 
the data in Figure HI to the spectral characteristics of H1426+428, and the Crab Nebula. 
The statistical correlation of the errors shown in this integral plot of excess events has been 
taken into account to derive optimal estimates of the spectral indices, differential and integral 
fluxes, and their errors. A summary of the results is given in Table [5J For the Crab Nebula a 
spectral index of 2.75 with an uncertainty of 6% (la) was derived. The peak of the differential 
detection rate of the Whipple telescope from a source with such a spectral index has been 
found to be around 360 GeV (2001 observing season). Based on the small Crab Nebula 
data-set, the spectral index as well as the differential and integral fluxes have been derived, 
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and are found to be consistent, at the 2<r level, with our previous observations (IHillas et al. 



19981 ). Analysis of the H1426+428 data-set indicates a substantially steeper spectral index 



of 3.55 with a 13% relative error. For this spectral index and the 366 digital counts size 
cut applied in the data analysis, the differential photon detection rate peaks at ~ 280 GeV. 
Table [5] shows estimates of the fluxes at 390 GeV assuming the same spectral index. This 
energy was chosen by increasing the size cut until the signal to noise ratio decreased to ~ 
4.5 cr; below this value the signal was not considered strong enough for reliable calculation 
of the H1426+428 flux. 
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Table 5. Results of the spectral analysis of H1426+428 and Crab Nebula data-sets. 



Parameter 


Crab 


H1426+428 


H1426+428 


x 2 


20.60 


22.72 


22.72 


d.f. 


16 


16 


16 


T(s) 


14.7 x 10 3 


117.3 x 10 3 


117.3 x 10 3 


a 


2.75 


3.55 


3.55 


E p (GeV) 


360 


280 


390 


F p (cm" 2 s _1 ) 


1.10 x 10~ 10 


2.04 x 10~ n 


8.76 x 10- 12 


dF p (cm" 2 s" 1 TeV- 1 ) 


5.34 x 10~ 10 


1.86 x 10~ 10 


5.73 x 10- 11 


vF v (erg cm -2 s -1 ) 


1.11 x 10~ 10 


2.33 x 10~ n 


1.39 x 10~ n 


Aa 

a 


0.06 


0.13 


0.13 


AF P 
F„ 


0.07 


0.17 


0.23 


AdFp 
dFp 


0.11 


0.25 


0.17 



Note. — a is the spectral index of the differential flux. E p denotes 
the energy at which the differential detection rate of photons from 
a given source, with given observation and data analysis conditions, 
peaks. F p and dF p are estimates of the integral and differential fluxes 
at E p . 
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Fig. 4. — Integral excess events observed by the Whipple telescope from the directions 
of H1426+428 (crosses) and the Crab Nebula (open circles) during 2001 as a function of 
integrated Cerenkov light in the shower image. Exposure on the Crab Nebula was adjusted 
to match the total excess of H1426+428 at the lowest size cut applied, 366 digital counts. 
One photoelectron corresponds to ~ 3.6 digital counts. 
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The observed energy fluxes of H1426+428 at 280 and 390 GeV are shown in Fi gure [5] 
together with the predicted SED for this source suggested by lCostamante et al.l (120011 ). This 
figure also shows the tentative detection (+3.1 a) made by the Whipple 10m telescope during 
2000. Operating at a substantially higher peak response energy at that time, the Whipple 
telescope observed this source at E pea k ~ 430 GeV. Based on these data an estimate of the 
energy flux of 0.86 ± 0.33 x 10~ n erg cm~ 2 s _1 has been derived. Due to the very weak 
signal however, we have not attempted to find the spectral index of H1426+428 for this 
data-set. As a result, the indicated error is only statistical and it does not include errors 
associated with the unknown spectral properties of the source. The 430 GeV peak energy 
itself is subject to a similar uncertainty. 
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Fig. 5.— SED of H1426+428 from ICostamante et all (boOll ). along with the observational 
results obtained with the Whipple telescope. The 280 GeV (log v = 25.83) and 390 GeV 
(log v = 25.97) data points, derived from 2001 data-set, suggest a differential spectral index 
of the energy flux equal to -3.55. The 430 GeV (log v = 25.99) point is a detection by the 
Whipple telescope during 2000. This result is subject to a large uncertainty associated with 
the weak detection and unknown spectral properties of the source. The EGRET upper limit 
above 100 MeV of 0.7 x 10~' cm -2 s _1 is als o shown as are the fluxes from various x-ray, 
optical, and radio experiments, all taken from ICostamante et al. J200lh . 
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Discussion 



Observations of H1426+428 with the Whipple 10m telescope since 1995 have returned 
consistently positive excesses although the statistical significance has not always been high. 
H1426+428 w as on the very short list of predicted TeV emitters published by the BeppoSAX 
collaboration (jCostamante et al.ll2000al ) and was also singled out as the most probable TeV 
emitter amongst these due to the high frequency of its synchrotron peak. It was observed as 
part of a BL Lac survey at the Whipple Observatory between 1995 and 1998, and was singled 
out for more extensive observations between 1999 and 2001. The strongest evidence for VHE 
emission was recorded during the 2000 and 2001 observing seasons, with signals at the 3.1a 
and 5.5a levels, respectively, being detected. Although it is difficult to combine results from 
different configurations of the camera, the combined significance of the results from 2000 and 
2001 is clearly equivalent to a greater than 5a detection, the usual level required to establish 
the existence of a new TeV source. 

There have been two occasions (in 1999 and in 2000) when there was marginal evidence 
for a transient signal from H1426+428. In 1999, after 2.3 hours of observations (over two 
nights) a signal at the +3. la level was recorded while in 2000, a signal at the +3.7a level 
was recorded after 2.2 hours of observations (on a single night). 

After 18.1 hours of observations on H1426+428, the CAT collaboration r eported a 3 a 



upper limit of 3.54 x 10 cm s at an energy threshold of 250 GeV (jPironl 120001 ) . 



These observations were carried out at the beginning of February 1999 in response to an 
announcement from the BeppoSXX and the RXTE/ASM teams, which revealed H1426+428 
to be transitioning to a high state. 

If independently verified, the observations reported here represent an important addition 
to the catalog of TeV-emitting BL Lacs (Table [6]). That the source was predicted to be a 
TeV 7-ray source based on its x-ray spectrum is important in that it signifies the maturity 
of the observational techniques and the theoretical understanding of BL Lacs. It reinforces 
the symbiosis between observations at x-ray wavelengths, particularly in hard x-rays, and 
those at TeV energies, particularly those with good sensitivity below energies of 1 TeV. 
The existence of a population of sources whose most prominent emission is at energies of 
10-100 keV and 300-1000 GeV points to a fruit ful overlap between the next generation of 
ground-based atmospheri c Cerenkov telescopes (IWeekes et al.ll200fl ) and the future hard x- 
ray experiment, EXIST (IGrindlay et al.ll2000l ). It is remarkable also in that this source, 
like Mrk 501 an d 1E2344+514, is not included in the 3rd EGRET Catalog, supporting the 
interpretation of iGhisellinil (119991 ) that in blazars sufficiently powerful to be detected readily 
by EGRET, the spectrum will terminate at lower energies and not extend to the TeV domain. 
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Table 6. Extra-galactic TeV Sources (jCatanese fc Weekeslll999l ) 



Source Type 



Discovery 



Group 



EGRET 



Markarian 421 HBL 



Markarian 501 HBL 



1ES 2344+514 HBL 



1ES 1959+650 HBL 



PKS 2155-304 HBL 



H1426+428 HBL 



3C66A 



LBL 



BL Lacertae LBL 



0.031 1992 Whipple yes 

fPunch et al. 19921 
0.034 1995 Whipple yes 

fQuinn et al. 19961 
0.044 1997 Whipple no 

fCatanese et al. 1998) 
0.048 1999 Telescope Array no 

fNishivama et al. 20001 
0.116 1999 Durham yes 

CChadwick et al. 19991 
0.129 2001 Whipple no 

(Horan et al. 2000, 2001a, 2001b) 
0.444 1998 Crimea yes 

fNeshpor et al. 19981 
0.069 2001 Crimea yes 

fNeshpor et al. 20011 
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Although the properties of H1426+428 reported here are scanty in comparison with the 
other, better studied TeV BL Lacs, they agree in principle with the characteristics now well 
established for Mrk 421 and Mrk 501 (and to a lesser extent, 1ES 2344+514). Although the 
flux from H1426+428 is low, there is some evidence for time variability There is also an 
indication of a hard 7-ray spectrum. There is no evidence for a correlation of x-ray and TeV 
7-ray fluxes (unlike Mrk 501) but the sensitivity of observations to date is limited and even 
in the stronger sources this correlation has been shown to be complicated. 



Ghisellinil (119991 ) has proposed that there is a continuous sequence of BL Lacs with the 
peak of the synchrotron spectrum increasing as the luminosity decreases. H1426+428 is 
another example of an "extreme" BL Lac characterized by a synchrotron spectrum that may 
peak near 100 keV and a relatively weak luminosity. These objects are the best candidates 
for TeV emission based on Compton-synchrotron models of BL Lac jets. However, the 
observed variations in sources such as Mrk 501 cannot be accounted for by simple one- zone 
homogeneous self-Compton models and may require another source of optical target photons. 
It will require the detection of a greate r population of sources such as H1426+428 to test 



these models fully. It has been noted (jGhisellinil Il998l ) that these TeV-emitting BL Lacs 



are also characterized by a relatively strong radio luminosity and that the radio luminosity 
may be a measure of the density of additional seed photons at infrared-optical wavelengths 
needed to explain the TeV emission. 

This is the most distant of the TeV-detected BL Lacs classified as HBL, and hence it 
has promising implications for the detection of more BL Lacs at z > 0.1. Stronger detections 
of such sources, which allow an accurate measure of the TeV energy spectrum, may place 
significant limits on the density of the intergalactic background light. The steep spectrum 
derived here is consistent with many models of intergalactic absorption but could also be 
intrinsic to the source. 

Interpretation of the theoretical prediction of the intrinsic spectrum of H1426+428 and 
the observational results reported in this paper are complicated by the possible strong at- 
tenuation of the high energy photons by the diffuse intergalactic infrared background. The 
appearance of such a cutoff in the spectra of TeV extr agalactic sources was suggested in on e 



of the pioneering works on TeV 7-ray absorption, by IStecker. de Jager fc Salamonl (119921 ). 
for the quasar 3C 279 (z = 0.54). A later paper (V. V. Vassiliev et al., in preparation) will 
discuss this subject in detail in relation to H1426+428; here we note only that at 280 GeV 
the attenuation optical depth can be anywhere between z of 0.1 and 1.0 depending on the 
extreme upper and lower limits known for the density of the extragalactic background light. 
Due to this potentially large absorption effect it is possible, for example, that the intrinsic 
spectral index of this source is softened by ~ 1.75 to produce the observed value of 3.55, if 



31 



the intrinsic properties of H1426+428 are similar to those of Mrk 501 (jVassiliev et al.ll200ll ). 
It remains certain however, that being the most distant HBL detected at sub-TeV energies, 
the observed spectral properties of H1426+428 and the understanding of its intrinsic char- 
acteristics from multiwavelength observations will provide the most constraining data for 
extragalactic background light studies. 
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A. Data Log for 2000 and 2001 Observations 
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Table 7: Data log for the 2000 and 2001 H1426+428 Observations 



MJD Start Time [Mode] 



51577.48031 


[p] 


51672 


28898 


[t] 


51941.43213 


[t] 


51970.44161 


[t] 


51994.37299 


[P] 

Ll J 


51577.52101 


[t] 


51687 


17778 


[t] 


51941.52991 


W 


51970.49679 


[t] 


51994.41703 


[P] 


51578.49144 


[p] 


51687 


19735 


[t] 


51942.43760 


W 


51970.51712 


[t] 


51995.49070 


[P] 


51579.47485 


[p] 


51687 


21689 


[t] 


51944.49363 


W 


51971.43684 


[t] 


51996.43112 


W 

L J 


51579.51770 


[t] 


51687 


23145 


[t] 


51957.37136 


[P] 


51971.45942 


[t] 


51996.47550 




51581.48048 


[p] 


51688 


17238 


[p] 


51957.39284 


W 


51972.44811 


[t] 


51997.35435 


[P] 

L-L J 


51585.50929 


[t] 


51689 


17450 


[p] 


51959.39676 


[P] 


51972.47175 


[t] 


51997.44241 


L J 


51586.50519 


[t] 


51690 


17385 


[p] 


51959.42716 


W 


51986.36833 


[p] 


51997.46204 


W 

L J 


51604.43009 


[t] 


51692 


20786 


[t] 


51959.44517 


W 


51988.33540 


[p] 


51998.40471 


[t] 


51605.40456 


[p] 


51692 


24942 


[p] 


51960.37838 


[P] 


51988.35239 


[t] 


51998.44546 


[P] 


51607.43234 


[p] 


51693 


16773 


[p] 


51960.40149 


W 


51989.38498 


[t] 


51999.41745 


[P] 


51617.42927 


[p] 


51694 


17954 


[t] 


51960.44430 


[P] 


51989.43090 


[p] 


51999.45334 


[P] 


51632.35804 


[t] 


51694 


19409 


[t] 


51960.46341 




51990.39131 


[t] 


52000.40087 


[t] 


51632.37754 


[t] 


51694 


22455 


[p] 


51960.48046 


W 


51990.43151 


[p] 


52000.44380 


[P] 


51635.31200 


[t] 


51694 


26542 


[p] 


51960.50067 


W 


51990.46683 


[t] 


52020.23130 


[P] 


51635.33273 


[t] 


51694 


30629 


[p] 


51961.37826 


[P] 


51991.40168 


[p] 


52022.57165 


[P] 


51636.35315 


[p] 


51695 


17682 


[p] 


51961.39750 


W 


51991.41978 


[t] 


52046.22484 


[P] 


51637.30195 


[p] 


51695 


21768 


[p] 


51961.43929 


[P] 


51991.43755 


[t] 


52047.31743 


L J 


51638.35877 


[p] 


51695 


25854 


[p] 


51961.47573 




51991.46142 


[t] 


52049.27725 


W 


51641.30764 


[p] 


51696 


26489 


[t] 


51961.49286 




51991.47977 


[t] 


52050.19936 


[P] 


51642.36119 


[t] 


51698 


17804 


[t] 


51962.39036 


[P] 


51992.38284 


[p] 


52050.23721 


[P] 


51642.37579 


[t] 


51699 


17421 


[p] 


51962.40577 




51992.43252 


[p] 


52051.19595 


[P] 


51644.36277 


[t] 


51700 


23519 


[p] 


51962.43053 


[P] 


51992.47643 


[t] 


52051.26034 


[P] 


51645.36182 


[p] 


51701 


25211 


[p] 


51963.38181 




52275.18178 


[p] 


52053.23280 


[P] 


51659.25434 


[p] 


51702 


28313 


[t] 


51963.42351 




52275.18443 


[t] 


52055.26242 


[P] 


51659.29540 


[t] 


51702 


30271 


[t] 


51963.46618 




52275.18320 


[t] 


52056.30267 


[P] 


51667.28374 


[p] 


51720 


18600 


[t] 


51963.48589 




52275.18199 


[t] 


52071.24343 


W 


51669.26859 


[p] 


51721 


18547 


[p] 


51965.45873 


[t] 


52275.18201 


[t] 


52072.21843 




51672.26933 


[t] 


51940 


53791 


[t] 


51967.46276 




52275.18178 


[t] 







Note. — The number in parentheses after the MJD gives the mode in which the data were taken: p - 

PAIRS, t - TRACKING. 
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